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■ ABSTRACT 

^ Aims. Our aim is to show how different hypotheses about Type la supernova progenitors can affect Galactic chemical evolution. Supernovae 
$H ■ la, in fact, are believed to be the main producers of Fe and the timescale with which Fe is restored into the interstellar medium depends on the 
[/3 ' assumed supernova progenitor model. This is a way of selecting the right progenitor model for supernovae la, a still debated issue. 

Metliods. We include different Type la SN progenitor models, identified by their distribution of time delays, in a very detailed chemical 
evolution model for the Milky Way which follows the evolution of several chemical species. We test the single degenerate and the double 
I degenerate models for supernova la progenitors, as well as other more empirical models based on differences in the time delay distributions. 
^ 1 Resuits. We find that assuming the single degenerate or the double degenerate scenario produces negligible differences in the predicted [O/Fe] 
' vs. [Fe/H] relation. On the other hand, assuming a percentage of prompt (exploding in the first 100 Myr) Type la supernovae of 50%, or that the 
f — maximum Type la rate is reached after 3-4 Gyr from the beginning of star formation, as suggested by several authors, produces more noticeable 

■ effects on the [O/Fe] trend. However, given the spread still existing in the observational data no model can be firmly excluded on the basis of 
' only the [O/Fe] ratios. On the other hand, when the predictions of the different models are compared with the G-dwarf metallicity distribution, 

1^-^ , the scenarios with very few prompt Type la supernovae can be excluded. 

Conciusions. Models including the single degenerate or double degenerate scenario with a percentage of 10-13% of prompt Type la supernovae 
^\ produce results in very good agreement with the observations. A fraction of prompt Type la supernovae larger than 30% worsens the agreement 
. with observations and the same occurs if no prompt Type la supernovae are allowed. In particular, two empirical models for the Type la SN 
' progenitors can be excluded: the one without prompt Type la supemovae and the one assuming delay time distribution going like oc r"^. We 
conclude that the typical timescale for the Fe enrichment in the Milky Way is around 1-1.5 Gyr and that Type la supernovae should appear 
already during the halo phase. 

Key words. Galaxy: evolution, Galaxy: abundances, Stars: supernovae: general - Galaxy: solar neighbourhood 
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1. Introduction 

The Type la supernova (SN) rate, which depends on the SN 
la progenitor model and the star formation history, is a fun- 
damental ingredient in models of galactic chemical evolution. 
The progenitor model can be described simply by a delay time 
distribution (DTD) function, which is the distribution of the 
explosion times. In the pioneering work of Greggio & Renzini 
(1983a, hereafter GR83) there was, for the first time, an expres- 
sion for the Type la SN rate in the scenario of the single degen- 
erate model. In this scenario, SNe Type la arise from the explo- 
sion of a C-O white dwarf in a close binary system where the 
companion is either a red giant or a main sequence star (Whelan 
& Iben, 1973; Munari & Renzini, 1992; Kenyon et al.l993; 
Hachisu et al., 1996; 1999; Han & Podsiadlowski, 2004). The 



DTD of this model arises from the above described scenario 
and predicts the explosion times of binary systems with the ap- 
propriate characteristics. An alternative model for progenitors 
of Type la SNe was proposed by Iben & Tutukov (1984). In 
this scenario, two C-O white dwarfs merge after loss of an- 
gular momentun due to gravitational wave emission, and ex- 
plode because the final object reaches the Chandrasekhar mass. 
Tornambe & Matteucci (1986) formulated a Type la SN rate in 
this scenario and applied it to galactic chemical evolution mod- 
els. In the following, Matteucci & Greggio (1986) tested the 
GR83 rate by means of a detailed chemical evolution model of 
the Milky Way and interpreted the [a/Fe] ratios versus [Fe/H] 
as due to the delay in the Fe production from Type la SNe, 
thus confirming previous suggestions (Tinsley, 1979; Greggio 
& Renzini 1983b). In recent years other authors (Dahlen et 
al. 2004; Strolger et al. 2004; Mannucci et al. 2005, 2006; 



* email to: matteucc@oats.inaf.it 



2 



Matteucci & al.: Different Supernova Progenitors 



Scannapieco & Bildsten, 2005; Sullivan et al. 2006; Aubourg 
et al. 2008; Pritchet et al. 2008; Totani et al. 2008) have pro- 
posed Type la SN rates based on DTDs derived empirically. 
One example is the rate of Mannucci et al. (2005;2006): they 
suggested that two populations of progenitors of Type la SNe 
are needed to explain the dependence of the rates on the colors 
of the parent galaxy. The presence of Type la SNe in old, red, 
quiescent galaxies is an indication that part of these SNe origi- 
nate from old stellar populations. On the contrary, the increase 
of the rate in blue galaxies (by a factor of 30 going from (B-K) 
~2 to (B-K) ~ 4.5) shows that part of the SN la progenitors 
is related to young stars and closely follows the evolution of 
the star formation rate (hereafter SFR). Mannucci, Delia Valle 
& Panagia (2006, hereafter MVP06), on the basis of the previ- 
ously described relation between the Type la SN rate and the 
color of the parent galaxies, their radio power as measured by 
Delia Valle et al. (2005), and cosmic age, concluded that there 
are two populations of progenitors of Type la SNe. They have 
suggested (but see Greggio, Renzini & Daddi, 2008) that the 
current observations can be accounted for only if about half of 
the SNe la (prompt SNe Id) explode within 10** years after the 
formation of their progenitors, while the rest explode during a 
wide period of time extending up to 10 Gyrs (tardy SNe la). 
Matteucci et al. (2006) appUed this formulation of the Type la 
SN rate to chemical evolution models of galaxies of different 
morphological type. They concluded that a fraction of 50% of 
prompt Type la SNe worsens the agreement with abundance 
data, especially in the Milky Way, and suggested that prompt 
Type la SNe should indeed exist but their fraction should not be 
larger than 30%. They also found a possible scenario for SN la 
progenitors in the MVP06 framework: a possible justification 
for this strongly bimodal DTD can be found in the framework 
of the single degenerate model for the progenitors of Type la 
SNe. In fact, such a DTD can be found if one assumes that 
the function describing the distribution of mass ratios inside 
the binary systems, is a multi-slope function. In particular, this 
choice means that in the range 5-8 M© are preferred the sys- 
tems where M\ ~ M2, whereas for lower mass progenitors are 
favored systems where Mi » M2. Another suggestion came 
from Dahlen et al. (2004) and Strolger et al.(2004, 2005), on 
the basis of the observational cosmic Type la SN rate. They 
suggested a DTD with no prompt Type la SNe with a maxi- 
mum occurring at 3-4Gyr. This DTD has never been tested in 
chemical evolution models. More recently, a direct determina- 
tion of the DTD function has been reported by Totani et al. 
(2008), on the basis of the faint variable objects detected in 
Subaru/XMM-Newton Deep Survey (SXDS). They concluded 
that the DTD function is inversely proportional to the delay 
time, i.e. the DTD can be well described by a featurless power 
law (DTD oc r", with « ~ 1). A similar suggestion came from 
Pritchet et al. (2008), who suggested a DTD oc r"^*^-^. 

The maximum at which the the Type la SN rate occurs is a 
very important parameter which affects the shape of the [a/Fe] 
relation in galaxies. This timescale depends on the fraction of 
prompt/tardy Type la SNe and on the specific history of star 
formation in galaxies, as shown by Matteucci & Recchi (2001), 
being shorter than in the solar neighbourhood for ellipticals and 
longer for irregular systems. 



In this paper, for the first time we show the effects of adopt- 
ing different DTDs in a model for the chemical evolution of 
the Milky Way: in particular, we will focuse on three relevant 
observational costraints: i) the [O/Fe] vs.[Fe/H] relation in the 
solar vicinity for which many accurate data are now available, 
ii) the Type la SN rate at the present time and iii) the G-dwarf 
metallicity distribution. To do that we we adopt five different 
DTDs, two of which are empirically derived. 

The paper is organized as follows: in Section 2 the formu- 
lation for the SNIa rates is presented. In Section 3 the galactic 
chemical evolution model is described and the results are dis- 
cussed, and in Section 4 some conclusions are drawn. 

2. Type la supernovae 

2.1. Progenitors 

We recall here the most common models for the progenitors of 
Type la SNe proposed insofar: 

- The merging of two C-O white dwarfs (WDs), due to grav- 
itational wave radiation, which reach the Chandrasekhar 
mass (~ 1 .4-Mq) and explode by C-deflagration (Iben and 
Tutukov 1984). This is known as double-degenerate (DD) 
scenario. In the original paper of Iben & Tutukov the pro- 
genitor masses are in the range 5-9Mq to ensure two WDs 
of at least ~ Q.IMq in order to reach the Chandrasekhar 
mass. In Greggio (2005, hereafter G05) the range of pro- 
genitors is 2.O-8.OM0, so that different combinations of WD 
masses are allowed, although always leading to exceed the 
Chandrasekhar mass. The clock for the explosion in this 
scenario is given by the lifetime of the secondary plus the 
gravitational time delay. 

- The C-deflagration of a Chandrasekhar mass C-O WD after 
accretion from a non-degenerate companion (Whelan and 
Iben 1973; Munari and Renzini 1992; Kenyon et al. 1993). 
This model is known as the single-degenerate (SD) one. 
The main problem with this scenario is the narrow range of 
permitted values of the mass accretion rate in order to ob- 
tain a stable accretion, instead of an unstable accretion with 
a consequent nova explosion and mass loss. In this case, 
in fact, the WD never achieves the Chandrasekhar mass. In 
particular, Nomoto, Thielemann & Yokoi (1984) found that 
a central carbon-deflagration of a WD results for a high ac- 
cretion rate (M > 4 • 10"^ Mq yr ') from the secondary to 
the primary star (the WD). They found that ~ (0.6-0.7) Mq 
of Fe plus traces of elements from C to Si are produced 
in the deflagration, well reproducing the observed spectra. 
The clock for the explosion here is given by the lifetime of 
the secondary star. 

- A sub-Chandrasekhar C-O WD exploding by He- 
detonation induced by accretion of He-rich material from 
a He star companion (Tornambe and Matteucci 1987; 
Limongi and Tornambe 1991). 

- A further model by Hachisu et al. (1996; 1999) is based 
on the classical scenario of Whelan and Iben (1973), but 
they find an important metallicity effect. When the accre- 
tion process begins, the primary star (WD) develops an op- 
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tically thick wind which helps in stabilizing the mass trans- 
fer process. When the metallicity is low {{FejH] < -1), 
the stellar wind is too weak and the explosion cannot oc- 
cur. The clock for the explosion here is also given by the 
lifetime of the secondary star. However, one has to take into 
account that Type la SN progenitor systems do not form be- 
fore the gas has attained a threshold metallicity. Therefore, 
the chemical enrichment from Type la SNe is delayed be- 
cause of this effect. 

2.2. The assumed DTDs 

We have tested the effects of five DTDs: 

- The Matteucci & Recchi (2001, hereafter MROl) one. 
which is practically the same as GR83, and is based on 
the single degenerate model for Type I SN progenitors (see 
Matteucci et al. 2006 for details). In this DTD the fraction 
of prompt Type la SNe is 13%. This formulation has proven 
one of the best to describe the chemical evolution of most 
galaxies. The very first SNe la to explode in this scenario 
are systems made of an 8-1-8 M© which takes ~ 30-35 Myr. 

- The DTD suggested by G05 for the DD scenario: in partic- 
ular, the DTD characterized by the exponent of the distri- 
bution function of final separations /3a = -0.9, by the expo- 
nent of the gravitational delay time distribution /3g - -0.75 
and by a maximum nuclear delay time (lifetime of the sec- 
ondary star) T„ ^ = 0.4 Gyr. The reason for chosing these 
parameters resides in the fact that they ensure a not too flat 
distribution of separations, which seems to be required to 
reproduce the specific Type la SN rate in galaxies as a func- 
tion of colors (see Greggio and Cappellaro 2009). In this 
DTD the fraction of prompt Type la SNe is 10%. The very 
first SNe la to explode in this scenario are systems made of 
an 8h-8 Mq which takes - 30 - 35 Myr plus the minimum 
gravitational time delay (1 Myr, G05). 

- The DTD of MVP06. Here the fraction of prompt Type la 
SNe is 50%. The very first systems explode also after 30-35 
Myr. 

- The DTD proposed by Strolger et al. (2004, hereafter S04) 
with a maximum at 3-4 Gyr. Here the very first systems 
explode after 2.5 • 10^ years. Therefore, there are no prompt 
Type la SNe. 

- The DTD proposed by Pritchet et al. (2008, PHS08), with 
DTD oc f where the fraction of prompt Type la SNe is 
~ 4%. The one proposed by Totani et al. (2008) goes like 

and we did not test it since these authors have already 
shown that it is very similar to the DTD of G05 for the DD 
scenario. 

In Figure 1 we show the different DTDs including those men- 
tioned above. As one can see, while the DTDs of the SD and 
DD scenarios are similar, those of MVP06 and S04 and PHS08 
are quite different. We did not take into account the Hachisu 
et al. (1996; 1999) model since its effects have been already 
discussed in Kobayashi et al. (1998) and MROl. We did not 
consider also the sub-Chandra model because the correpond- 
ing SN rate was already computed by Tornambe and Matteucci 
(1987). Finally, we show the DTD relative to the DD model 
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Fig. 1. Various DTD functions normalized to their own maxi- 
mum value: the continuous blue line is the DTD of MROl; the 
long dashed green line is the DTD of G05 for the DD wide 
channel; the dotted magenta line is the DTD for the DD close 
channel of G05; the dashed red Une is the DTD of MVP06; the 
short dashed-dotted black line is the DTD of S04, the cyan long 
dashed-dotted hne is the DTD of PHS08. 

(close channel) but we will not use it for computing galactic 
chemical evolution, and the reason for this is that it is very sim- 
ilar to the DD (wide channel) DTD. 

2.3. The calculation of the Type la SN rate 

Once a DTD has been assumed, the Type la SN rate is com- 
puted according to Greggio (2005): 



Rla(t) 



A(t - T)iff(t - T)DTD{T)dT 



(1) 



where i//(t) is the SFR and A{t - r) is the fraction of binary 
systems which give rise to Type la SNe and in principle it can 
vary in time. Here we will assume A to be a constant. It is worth 
noting that the fraction A represents the fraction of binary sys- 
tems with those particular characteristics to give rise to Type 
la SNe relative to the whole range of star masses (O.l-IOOM©). 
The time t is the delay time defined in the range (t,, Tx) so that: 



DTD(T)dT = I 



(2) 



where t, is the minimum delay time for the occurrence of 
Type la SNe, in other words the time at which the first SNe 
la start occurring, and is the maximum delay time. Clearly 
these delay times vary according to the assumed progenitor 
scenario: in particular, in the SD scenario is just the max- 
imum nuclear delay time, namely the lifetime of the smallest 
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secondary star, whereas in the DD model is the sum of the max- 
imum nuclear delay time and the maximum gravitational delay 
time. 

Finally, ka is the number of stars per unit mass in a stellar 
generation and contains the IMF. 
In particular: 

ka = I (pim)dm (3) 

•JniL 

with the normalization condition for the IMF being: 

rmu 

I m4>im)dm = 1 (4) 

where nii - 0.1 Mq and mu = IOOMq and define the whole 
range of existence of the stars. 

3. The model for the Milky Way 

The model we adopt for computing the chemical evolution of 
the Milky Way is based on the original two-infall model by 
Chiappini et al (1997; 2001). In particular, they assumed that 
the first infall episode was responsible for the formation of the 
halo and thick-disk stars that originated from a fast dissipative 
collapse. The second infall episode formed the thin-disk com- 
ponent, with a timescale much longer than that of the thick-disk 
formation. The timescale for the formation of the halo-thick 
disk is 1-2 Gyr, while the timescale for the thin disk is 7 Gyr 
in the solar vicinity. The authors included in the model also 
a threshold in the gas density, below which the star formation 
process stops. The existence of such a threshold value is sug- 
gested by observations relative to the star formation in external 
disk galaxies (Kennicutt 1998, but see Boissier et al. 2007). The 
physical reason for a threshold in the star formation is related 
to the gravitational stability, according to which, below a criti- 
cal density, the gas is stable against density condensations and, 
consequently, the star formation is suppressed. 
In the two-infall model the halo-thick disk and the thin disk 
evolutions occur at different rates and they are independent, 
mostly as a result of difi'erent accretion rates. With these pre- 
cise prescriptions it is possible to reproduce the majority of the 
observed properties of the Milky Way and this shows how im- 
portant is the choice of the accretion law for the gas coupled 
with the SFR in the Galaxy evolution. Some of the most im- 
portant observational constraints are represented by the various 
relations between the abundances of metals (C, N, a-elements, 
iron peak elements) as functions of the [Fe/H] abundance (the 
most common tracer of metallicity) and by the G-dwarf metal- 
Ucity distribution. The Galactic disk is approximated by a se- 
ries of concentric annuli, 2 kpc wide, without exchange of mat- 
ter between them. In this model the thin disk forms "inside- 
out", in the sense that the timescale for the disk formation in- 
creases with galactocentric distance. This choice was dictated 
by the necessity of reproducing the abundance gradients along 
the Galactic disk (see Chiappini et al. 2001). 

The time-evolution of the fractional mass of the element i in 
the gas within a galaxy, G„ is described by the basic equation: 

Gi = -m^iit) + Ri(t) + (Gdinf - (Gdout (5) 



where G,(?) = Mg{t)Xi{t)IM,ot is the gas mass in the form 
of an element ; normalized to a total fixed mass M,oi and 
G{t) = Mg(t)/M,ot is the total fractional mass of gas present 
in the galaxy at the time t. For modelling spiral disks the 
masses are substituted by the surface mass densities. The quan- 
tity Xi(t) = Gi{t)IG{t) represents the abundance by mass of an 
element i, with the summation over all elements in the gas mix- 
ture being equal to unity. 

Ri{t) represents the returned fraction of matter in the form of an 
element i that the stars eject into the ISM through stellar winds 
and supernova explosions; this term contains all the prescrip- 
tions concerning the stellar yields and the supernova progenitor 
models. 

The assumed SFR {ij/it)) is a Schmidt (1955) law with a de- 
pendence on the surface gas density {k = 1.5, see Kennicutt 
1998) and also on the total surface mass density (see Dopita & 
Ryder 1994). In particular, the SFR is based on the law origi- 
nally suggested by Talbot & Arnett (1975) and then adapted by 
Chiosi(1980): 

/2(r,f)2„a,(r,/)\'*"'' t 

where t is the time, r is the galactocentric distance and the 
constant v is the efficiency of the star formation process and is 
expressed in Gyr ': in particular, v - 2 G>t ' for the halo and 
1 Gyr ' for the disk (/ > 1 Gyr). The total surface mass den- 
sity is represented by S(r, t), whereas 2(ro, t) is the total sur- 
face mass density at the solar position, assumed to be r© = 8 
kpc (Reid 1993). The quantity ^gasir, t) represents the surface 
gas density. These choices of values for the parameters allow 
the model to fit very well the observational constraints, in par- 
ticular in the solar vicinity. A threshold gas density for the star 
formation in the thin disk of IM^pc'^ is adopted in all the mod- 
els presented here. It is worth noting that the SFR described by 
eq.(6) is equivalent to a simple Kennicutt' s law {SFR = vS*^^) 
for what concerns the evolution of chemical abundances, as 
shown in Colavitti et al. (2009). We have used here the for- 
mulation of eq.(6) for continuity with all our previous papers. 
The IMF is that of Scalo (1986) normalized over a mass range 
of 0.1-100 Mq and it is assumed to be constant in space and 
time. 

The nucleosynthesis prescriptions are common to all mod- 
els and are taken from: Woosley & Weaver (1995) yields for 
massive stars (those of Fe relative to the solar chemical compo- 
sition and those of O metallicity dependent since this is the best 
combination for reproducing the observations, as shown by 
Frangois et al. 2004); van den Hoeck & Groenewegen (1997) 
metallicity dependent yields for low and intermediate mass 
stars (0.8 < M/Mq < 8) and Nomoto et al.'s yields (1997) for 
Type la SNe. The choice of these particular yields resides in the 
fact that they give the best agreement with the observed abun- 
dances in the solar vicinity, as shown by Francois et al. (2004). 
Moreover, the O yields are those best known, since the major- 
ity of nucleosynthesis studies agree on this particular element. 
The yields of Fe from massive stars are still more uncertain 
and they vary from author to author: however, most of Fe is 
assumed to be produced by Type la SNe and it is known from 
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observational studies, which suggest that these SNe produce, 
on average, 0.5Mq of Fe (see Blanc & Greggio, 2008 and ref- 
erences therein). The assumed stellar lifetimes for all galaxies 
are those suggested by Padovani & Matteucci (1993). 

The two terms (Gi)i„f and (Gdout account for the infall of 
external gas and for galactic winds, respectively. For the Milky 
Way we assume no wind, therefore (G,)oh/ - 0.. 
The infall instead is an important physical process in the 
Galaxy. We assume the two-infall law of Chiappini et al. 
(1997): 

A(r, t) = a{r)e-'l''"^'"> + b{r)e-^'-'~"'>'">^''^ (7) 

where a{r) and b{r) are two parameters fixed by repro- 
ducing the total present time surface mass density along the 
Galactic disk. It is worth noting that this particular infall law 
has proven to be consistent with an accretion law derived in 
a cosmological context, by assuming that gas accretion follows 
the accretion of the dark matter halo for a galaxy like the Milky 
Way (Colavitti et al. 2008). Colavitti et al. (2008) tested various 
infall laws including cosmologically derived ones and costant 
and linear laws. These different infall laws influence the SFR 
which depends on the amount of gas present in the galaxy at 
any time and therefore the chemical abundances and the stel- 
lar metallicity distribution. They concluded that the best law 
in order to fit the main properties of the Milky Way is still eq. 
(7) which has also been adopted in other chemical evolution 
studies (e.g. Chang et al. 1999; Alibes et al. 2001). 

In the solar vicinity the total surface mass density cr,„,= 51 
+ 6 Mq pc"^ (see Boissier & Prantzos 1999). f,„aj: = 1.0 Gyr is 
the time for the maximum infall on the thin disk, th = 2.0 Gyr 
is the time scale for the formation of the halo thick-disk and 
To(r) is the timescale for the formation of the thin disk and it is 
a function of the galactocentric distance (formation inside-out, 
Matteucci and Frangois 1989; Chiappini et al. 2001). 

In Figure 2 we show the star formation rate as predicted 
by eq. (6) in the framework of the two-infall model. Note the 
gap in the star formation rate at about IGyr; this corresponds 
to the transition phase between the halo-thick disk and the thin 
disk phase. The gap is due to the fact that the disk is formed 
by means of the second infall episode and star formation can 
occur only if the gas density is above the threshold. Therefore, 
the gap occurs as a natural consequence of the two infall epi- 
sides and a threshold gas density for star formation. It is worth 
noting that our chemical results would not be much affected if 
we had considered only one infall episode and no threshold in 
the gas density. In this case, the SFR would be simply expo- 
nential without oscillations and there would be no gap between 
the thick and thin disk. The [O/Fe] vs. [Fe/H] (see next sec- 
tion) would differ only for the lack of a loop in the transition 
phase between thick and thin disk (see for example Matteucci 
& Francois 1989; 1992), as it is instead predicted by the two- 
infall model. 

4. Predicted [O/Fe] vs. [Fe/H] with differents DTDs 

In Figure 3 we show the predicted Type la SN rates for the so- 
lar vicinity in the framework of the two-infall model previously 
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Fig. 2. The star formation rate as predicted by the two-infall 
model in the solar neighbourhood: the first peak represents the 
star formation rate during the halo-thick disk phase. Note that 
there is a short gap in the star formation between the end of the 
thick disk phase and the beginning of the thin disk phase: this is 
due to the adoption of a threshold density for the star formation. 
Note that the oscillating behaviour in the last few Gyr is due to 
the fact that the threshold density is easily reached al late times 
and star formation stops for a while until the gas restored by 
dying stars increases again the interstellar gas density. 

described and for the considered DTDs. We have chosen the pa- 
rameter A in eq. (1) in order to reproduce the present time Typel 
a SN rate in the Galaxy (~ 0.3 S Ne century^^, Cappellaro et 
al. 1999). It is worth noting that a value of A~ 0.0025 for 
all the DTDs was chosen, in agreement with Matteucci et al. 
(2006), with the exception of the PHS08 DTD for which the 
parameter A=0.0002. The PHS08 DTD is, in fact, predicting 
the largest fraction of SNe la exploding with large delays (see 
Fig. 1). That means that a larger fraction of long-living progeni- 
tors (born several Gyr ago, when the SFR was much higher, see 
Fig. 2), relative to the other models, is exploding at the present 
time. Therefore, a smaller value of A is required to reproduce 
the present-time SN la rate in the Galaxy. The oscillating be- 
haviour shown by some Type la SN rates in the last couple of 
Gyr is clearly related to the oscillations in the star formation 
rate. In particular, models with a larger percentage of prompt 
Type la SNe show the oscillations at late times (SD model and 
MVP06 model), whereas the models of S04, G05 and PHS08, 
having null, 10% and 4% fractions of prompt Type la SNe, re- 
spectively, show negligible oscillations. 

As one can see, it is evident the effect of the gap in the star 
formation rate, as predicted by the two-infall model, for the 
DTDs of MVP06, SD, DD and PHS08 models. For the DTD 
of S04, there is no gap since the minimum delay in the occur- 
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Fig. 3. The predicted Type la SN rates in the solar vicinity in 
the framework of the two-infall model. The five curves refer 
to different DTDs. In particular, the dotted line represents the 
DTD of MVP06; the continuous line represents the SD model 
with the DTD from MROl and the short dashed line the DD 
model with the DTD of G05 for the DD wide channel. Finally, 
the long dashed line represents the DTD suggested by S04 and 
the dashed-dotted line the DTD of PHS08. The rates are ex- 
pressed in SN per century, and are all normalized to reproduce 
the present time Type la SN rate in the MiUcy Way. 



rence of Type la SNe is much longer than for the other DTDs. 
There are no prompt Type la SNe and most of the SNe la start 
occurring only after the gap in the star formation. In Figure 4, 
5 and 6 we show the effect of different DTDs on the [O/Fe] 
vs. [Fe/H] diagram. As it is well known this diagram is inter- 
preted on the basis of the "time-delay model". The time delay 
model suggests that the [O/Fe] in the halo stars ([Fe/H]< -1.0) 
is slowly varying with [Fe/H] and reflects mainly the behaviour 
of the O/Fe ratio in massive stars: therefore, one expects a slight 
variation with the stellar mass. For [Fe/H]> -1.0, the [O/Fe] ra- 
tio starts declining as a consequence of the increase of the Fe 
production from Type la SNe. Clearly, the DTD is fundamen- 
tal in determining the moment, and therefore the [Fe/H] value, 
at which this process becomes important. We remind that this 
effect does not depend on details of the star formation history 
but mainly on the time delay between the chemical enrichment 
from SNe II and SNe la. On the basis of what just said, it is 
therefore important to see which DTD best fits the abundance 
data. In Figure 4 we show the results for the DTDs related to 
the SD and DD (wide channel) models. As one can see, the dif- 
ferences between the results obtained by means of the SD and 
DD models are small, as it was expected by the comparison of 
the DTDs connected to these scenarios (see also Valiante et al. 
2009). In particular, the amount of prompt Type la SNe in the 



Fig. 4. The predicted [O/Fe] vs. [Fe/H] in the solar vicinity 
for two different DTDs: dashed line (DTD of SD model from 
MROl), dashed-dotted (DTD of G05 for the wide DD channel). 
The model results are normalized to their own predicted solar 
values. The data are taken from the compilation of Frangois et 
al. (2004). The solar position is marked by a red dot. 



SD and DD scenarios is very similar, namely 13% and 10%, 
respectively, as opposed to the DTD of MVP06 where the per- 
centage of prompt SNe la is 50% . This difference is evident in 
the [O/Fe] ratios predicted by the MVP06 DTD (Fig.5) which 
are declining before and faster than for the SD and DD models. 
In fact, in the model with MVP06 DTD the [O/Fe] knee occurs 
at [Fe/H] — 1.8 dex, whereas for the results obtained with the 
SD and DD DTDs the knee occurs at [Fe/H] ~ -0.8 dex (see 
Table 1). 

In Figure 5 we show the predicted [O/Fe] trend for the DTD 
suggested by S04 where the delay in the appearence of the first 
Type la SNe is much longer than in the previous cases. As one 
can see from the figure, the [O/Fe] ratio in this case shows a 
longer plateau than in the other models, and the knee in the 
[O/Fe] ratio occurs at [Fe/H]~ -0.6 dex. It is worth noting that 
the cosmic epoch at which the knee in the [O/Fe] ratio occurs 
in the models with the SD and DD rates is ~1 Gyr, whereas in 
the case of Strolger et al. DTD this time is 4-5 Gyr. In Figure 5 
we also show a comparison between the results obtained with 
a DTD oc f-o 5 (PHS08) and the DTD of MVP06. As one can 
see, the PHS08 DTD produces results not very different from 
the SD model but it predicts a too low solar Fe abundance 
([Fe///]Q=-0. 16dex) and a too high solar [0/Fe]Q=Q.22dex, as 
opposed to all the other models where the solar [Fe/H]^ and 
[0/Fe]Q are very close to zero (see Table 1 where we show the 
predicted solar ratios). Note that the model results are always 
normalized to their own predicted solar values. 
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Table 1. Predicted [O/Fe]^ and [Fe/H]^ values at the time of formation of the solar system plus the value of [FelH]i^nee at 
which the knee for each model occurs. Normalization to the solar values of Asplund et al. (2005). 



Model 






[FelH]knee 


SD 


0.095 


-0.013 


-0.8 


DD 


0.095 


-0.038 


-0.8 


MVP06 


0.07 


-0.003 


-1.8 


S04 


0.06 


0.01 


-0.6 


PHS08 


0.22 


-0.16 


-0.45 




Fig. 5. In this figure we show the predicted [O/Fe] when the 
Strolger et al. DTD is assumed (dotted line) compared with 
the predictions from the models with MVP06 DTS (continuous 
Une) and the model with the DTDoc f " ^ (PHS08) (long dashed 
line. The model results are normalized to their own predicted 
solar values. The data are the same as in Fig. 4. 



Finally, in Figure 6 we show a zoomed comparison of all 
the studied DTDs in the [Fe/H] range -1.5-0.5 dex: here the 
different effects produced by the different assumptions about 
Type la SN progenitors are more evident. A larger number of 
prompt Type la SNe clearly produces lower [O/Fe] ratios in 
this [Fe/H] range with respect to models depressing the num- 
ber of prompt SNe in favor of tardy ones. Although the spread 
existing in the data prevents us from drawing firm conclusions, 
from this comparison it turns out that the S04 model tends to 
predict too high values of [O/Fe] in the range -1.5 < [Fe/H]< 
-0.2.On the other hand, the model with MVP06 DTD predicts 
the lowest [O/Fe] ratios in the range -1.8 < [Fe/H]< -1.0. From 
numerical experiments, we can say that the MVP06 DTD with 
~ 30% of prompt SNe would fit better the data, and it would be 
still compatible with the findings of MVP06. It is also interest- 
ing to note that the DTDs influence the [Fe/H] at which the gap 



Fig. 6. Comparison of all the models with all the studied DTDs: 
the short dashed curve refers to the DTD of MROl; the contin- 
uous curve refers tothe DTD of MVP06; the dotted line rep- 
resents S04 DTD; the long dashed curve refers to the DTD of 
PHS08; the dashed-dotted curve refers to the DTD of the DD 
scenario (wide channel) of G05. The model results are normal- 
ized to their own predicted solar values. The data are the same 
as in Fig.4. 



in the star formation between the halo-thick disk and the thin 
disk phases is visible. This gap produces a loop in the [O/Fe] 
behaviour and is naturally produced by assuming a threshold 
gas density in the star formation. Observationally, a gap in the 
abundance ratios has been claimed by Fuhrmann (1998) and 
Gratton et al. (2000). In particular, the gap was observed at a 
[Fe/H] between -1.0 and -0.5 dex, corresponding to the tran- 
sition phase between the thick and the thin disk. If confirmed, 
this gap can impose constraints on the history of star formation. 
In our plots the gap occurs at larger [Fe/H] values going from 
the SD to PHS08, DD and MVP06 DTD (see Figure 6). For 
the S04 DTD, the gap is not visible and the reason is that Type 
la SNe start exploding in a non negligible number only after 
the gap. Therefore, the effect of the halt of oxygen production. 
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-1 -0.5 0.5 

[Fe/H] 



Fig. 7. The G-dwarf metallicity distribution in the solar vicin- 
ity. The data are from the Geneva-Copenhagen Survey of the 
Solar Neighbourhood, Nordstrom et al. (2004) (thick line his- 
togram) and the data by J0rgensen (2000) (thin-line histogram). 
The models refer to: S04 DTD (dotted Une), PSH08 DTD (long 
dashed Hne), DD model (dashed-dotted Hne), MVP DTD (con- 
tinuous line) and SD model (short dashed Une). 



while the Fe is still produced by the Type la SNe exploding dur- 
ing the gap, is lost. This effect produces in the other models the 
loop in the [O/Fe] ratio. If the gap really exists this would fur- 
ther strenghten the conclusion that prompt Type la SNe should 
exist and be a non negligible fraction of all SNe la. 
In summary, the SD and DD models together with the bimodal 
distribution of MVP06 all produce acceptable results, although 
a 50% of prompt Type la SNe appears too large. The DTD from 
PHS08 and S04 are probably less likely although the spread in 
the data prevents us from drawing firm conclusions. 
Therefore, we checked our results on a very important observa- 
tional constraint: the G-dwarf metallicity distribution. In Figure 
7 we show the predictions for the G-dwarf metallicity distribu- 
tion of all our models compared with observational data. This 
Figure shows that two models can be ruled out: the PHS08 
and the S04. Both of these models contain empirically derived 
DTDs. They predict too many metal poor stars below [Fe/H]< 
-1.0 dex (the so-called G-dwarf problem), and this is due to the 
fact that in both cases the Fe abundance increases more slowly 
that in the other three models (SD, DD, MVP06) thus produc- 
ing too many metal poor stars. The reason for that resides in 
the small number of prompt Type I a SNe in the PHS08 and 
S04 models. In particular, the S04 model shows two distinct 
peaks in the G-dwarf metallicity distribution and this is due to 
the fact that at the beginning of the thin disk formation only 
SNell enrich the gas in Fe and the bulk of Fe comes with a 
large delay, thus creating the second peak. On the other hand. 



the other models all produce acceptable results, with the DD 
and MVP06 models being the best. Note that all the models 
predict less metal rich stars than observed by Norstrom et al. 
(2004) which is the largest and more accurate survey of G and 
F dwarfs. This is a problem for pure models of chemical evo- 
lution which do not take into account the possibility of stel- 
lar migration. In fact, because of interaction between stars and 
transient spiral density waves (Roskar et al. 2008), stars born in 
the inner part of the disk can be scattered at larger galactocen- 
tric distances and this could be the explanation for the existence 
of metal rich stars in the solar vicinity. 

5. Conclusions 

We have studied the effects of different DTDs on the [O/Fe] 
vs. [Fe/H] diagram, where the effects of the time delay in the Fe 
production by Type la SNe are particularly important, and on 
the G-dwarf -[Fe/H] distribution. We considered the DTDs re- 
lated to the two main progenitor models for SNe la, namely the 
SD and DD models. Then we considered other DTDs derived 
empirically from the study of SN rates in galaxies (MVP06, 
PHS08) and the cosmic SN rate (S04), which are not di- 
rectly related to particular progenitor models. We computed the 
chemical evolution of the solar vicinity by adopting a success- 
ful model of Galactic chemical evolution, the so-called two- 
infall model. Our main results can be summarized as follows: 

- Delay time distributions containing a not negligible number 
of prompt Type a SNe (namely those SNe exploding inside 
100 Myr since the beginning of star formation) are neces- 
sary to best fit the [O/Fe] ratio. This indicates that the first 
Type la SNe occurring in galaxies can explode as soon as 
after 35-40 Myr since the beginning of star formation. This 
coiTesponds to the lifetime of a 8M0 star which is also the 
maximum mass for the progenitors of C-O white dwarfs. As 
a consequence of this, the very first Type la SNe already ap- 
pear during the halo phase. 

- The bimodal DTD of MVP06 predicts a very high fraction 
of prompt Type la SNe (~ 50%). This fraction is probably 
too high since it depresses the [O/Fe] ratio in the [Fe/H] 
range -1.8-0.8 dex and produces a knee in the [O/Fe] ratio 
at a too low metallicy (the observed knee is at [Fe/H] ~ 
-1.0 dex). A smaller fraction of promt SNe la would be 
acceptable (see also Matteucci et al. 2006). 

- The DTD of S04, which predicts the maximum of Type la 
SNe occurring at 3.4 Gyr since the beginning of star for- 
mation, overproduces the [O/Fe] ratio in the [Fe/H] range 
-1.0- -0.2 dex and predicts a knee in the [O/Fe] ratio oc- 
curring at a too high [Fe/H] ~ -0.6 dex. Although the spread 
in the data prevents us from drawing firm conclusions the 
best DTDs from the point of view of the [O/Fe] diagram 
are the one related to the SD and DD cases. 

- The G-dwarf metallicity distributions predicted by the 
models with different DTDs give a more clear indication 
on the DTDs which should be discarded and those are the 
two empirical ones derived by S04 and PSH08. They pro- 
duce too many disk stars with metallicity [Fe/H]< -1.0 dex. 
Therefore, we conclude that the DTDs related to the two 
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classical models for Type la SN progenitors (SD and DD), 
and in particular that related to the DD model (wide chan- 
nel,) are the best to reproduce the observational data. Since 
these two models include prompt Type la SNe even if in 
a smaller proportion than the bimodal DTD of MVP06, 
we conclude again that such prompt SNe are necessary 
to reproduce the chemical evolution of the solar vicinity. 
Probably a mixed scenario (SD+DD) would be the most 
reaUstic one (see Greggio et al. 2008). 
- On the basis of the results discussed above we can identify 
the typical timescale for Fe enrichment by Type la SNe in 
the solar vicinity, identified by the time at which the knee 
in the [O/Fe] ratio occurs, as 1-1.5 Gyr. 

6. Acknowledgments 

This work was partially supported by the Italian Space Agency 
through contract ASI-INAF 1/016/07/0. P.M., E.S. and S.R. 
acknowledge financial support from PRIN2007-MUR (Italian 
Ministry of University and Research), Prot.2007JJC53X_001. 
We also thanks an anonymous referee for his/her useful com- 
ments.S.R. acknowledges financial support from the FWF 
through the Lise Meitner grant M1079-N16. 

References 

AUbes, A., Labay, J.& Canal, R, 2001, A&A, 370, 1103. 
Asplund M., Grevesse N. & Sauval A. J., 2005, ASPC, 336, 
25 

Aubourg, E., Tojeiro, R., Jimenez, R., Heavens, A., Strauss, 

M. A., Spergel, D. N.,2008, A&A, 492, 631 

Blanc, G. & Greggio, L., 2008, NewA, 13 606 

Boissier S., Gil de Paz A., Boselli A., Madore B. P., Buat V., 

Cortese L., Burgarella D., Munoz Mateos J. C. et al., 2007, 

ApJSuppl., 173, 524 

Boissier, S., Prantzos, N. 1999, MNRAS, 307, 857 
Cappellaro, E., Evans, R., Turatto, M., 1999, A&A, 351, 459 
Chang, R. X., Hou, J. L., Shu, C. G.& Fu, C. Q., 1999, 

A&A,350, 38 

Chiappini C, Matteucci R & Gratton R., 1997, ApJ, 477, 765 
Chiappini C, Matteucci F. & Romano D., 2001, ApJ, 554, 
1044 

Chiosi C, 1980, A&A, 83, 206 

Colavitti E., Matteucci F. & Murante G., 2008, A&A, 483, 
401 

Colavitti E., Cescutti, G., Matteucci F. & Murante G., 2009, 

A&A, in press (astro-ph/081 1.3505) 

Dahlen, T., Strolger, L.G., Riess, A.G. et al. 2004, ApJ 613, 
189 

Delia Valle, M., Panagia, N., Padovani, P., Cappellaro, E., 

Mannucci, P, Turatto, M., 2005, ApJ 629, 750 

Dopita, M.A., Ryder,S.D,.1994, ApJ, 430,163 

Francois P., Matteucci P., Cayrel R., Spite M., Spite F. & 

Chiappini C, 2004, A&A, 421, 613 

Fuhrmann, K., 1998, A&A, 338, 161 

Gratton, R. G.; Carretta, E.; Matteucci, F.; Sneden, C. 2000, 

A&A, 358, 671 

Greggio, L., 2005, A&A 441, 1055 (G05) 



Greggio, L., Cappellaro, E., 2009, astro-ph/0902.0700 
Greggio, L., Renzini, A., 1983a, A&A, 118, 217 (GR83) 
Greggio, L., Renzini, A., 1983b, Mem. Salt, Vol. 54, p.311 
Greggio, L.,Renzini, A., Daddi, E., 2008, MNRAS, 388, 829 
Hachisu, 1., Kato, M. & Nomoto, K., 1996, ApJ 470, 97 
Hachisu, I., Kato, M. & Nomoto, K.,1999, ApJ 522, 487 
Han Z. & Podsiadlowski R, 2004, MNRAS, 350, 1301 
J0rgensen, B. R., 2000, A&A, 363, 947 
KennicuttR.C. Jr., 1998, ARA&A 36, 189 
Kenyon, S.J., Livio, M., Mikolajewska, J., Tout, C. A., 1993, 
ApJ, 407, L81 

Kobayashi, C, Tsujimoto, T, Nomoto, K., Hachisu, I. Kato, 
M., 1998, ApJ, 503, L155 

Iben, l.Jr. & Tutukov, A. 1984, ApJSuppl., 54, 335 
Limongi, M. Tornambe, A., 1991, ApJ, 371, 317 
Mannucci, F., Delia Valle, M., Panagia, N., 2006, MNRAS 

370, 773 (MVP06) 

Mannucci, P., Delia Valle, M., Panagia, N., Cappellaro, E., 
Cresci, G., Maiolino, R., Petrosian, A. & Turatto, M., A&A 

2005, A&A 433, 807 

Matteucci P. & Prangois R,1992, A&A, 262, LI 
Matteucci F. & Frangois R, 1989, MNRAS, 239, 885 
Matteucci P & Greggio, L., 1986, A&A, 154, 279 
Matteucci, P., Panagia, N., Pipino, A., Mannucci, P., Recchi, 
S. & Delia Valle, M., 2006, MNRAS 372, 265 
Matteucci, P& Recchi, S., 2001, ApJ, 558, 351 (MROl) 
Munari, U., Renzini, A., 1992, ApJ, 397, L87 
Nomoto, K, Thielemann, F.-K., Yokoi, K., 1984, ApJ, 286, 
644 

Nomoto, K., Iwamoto, K., Nakasato, N., Thielemann, P. K. et 
al. 1997, NuPhA 621, 467 

Nordstrom, B., Mayor, M., Andersen, J., Holmberg, J., Pont, 
P., Jrgensen, B. R., Olsen, E. H., Udry, S., Mowlavi, N., 2004, 
A&A, 418, 989 

Padovani, P, Matteucci, P, 1993, ApJ, 416, 26 

Pritchet, C. J., Howell, D. A., Sullivan, M., 2008, ApJ, 683, 

L25 (PHS08) 

Reid M. J., 1993, ARA&A, 31, 345 

Roskar, R., Debattista, V. P., Quinn, T. R., Stinson, G.S., 

Wadsley, J., 2008, ApJ, 684, L79 

Scalo J. M., 1986, PCPh, 11, 1 

Scannapieco, P., Bildsten, L.,2005 ApJ 629, L85 

Schmidt, M., 1955, ApJ, 129, 243 

Strolger, L.G., Riess, A.G., Dahlen, T. et al., 2004, ApJ 613, 
200 (S04) 

Strolger, L.G., Riess, A.G., Dahlen, T. et al., 2005, ApJ 635, 
1370 

Sullivan, M., Le Borgne, D., Pritchet, C. J., Hodsman, A., 

Neill, J. D., Howell, D. A., Carlberg, R. G., Astier, P., 

Aubourg, E., Balam, D. & al., 2006, ApJ, 648, 868 

Talbot R. J. Jr. & Arnett W. D., 1975, ApJ, 197, 551 

Tinsley, B.M., 1979, ApJ, 229, 1046 

Tornambe, A., Matteucci, P,1986, MNRAS, 223, 69 

Tornambe, A., Matteucci, P,1987, ApJ, 318, L25 

Totani, T. Morokuma, T, Oda, T., Doi, M., Yasuda, N. 2008, 

PAS J, 60, 1327 

Valiante, R., Matteucci, P. Recchi, S., Calura, P., 2009, New 
Astronomy in press (astro-ph/0807.2354v2) 



10 



Matteucci & al.: Different Supernova Progenitors 



van den Hoek L. B. & Groenewegen M. A. T., 1997, A&AS, 

123, 305 

Whelan J. & Iben I. Jr., 1973, ApJ, 186, 1007 

Woosley S. E. & Weaver T. A., 1995, ApJSuppl. 101, 181 



